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ON OPTIMAL ARBITRAGE 

By Daniel Fernholz and Ioannis Karatzas 1 

University of Texas at Austin and INTECH 

In a Markovian model for a financial market, we characterize 
the best arbitrage with respect to the market portfolio that can be 
achieved using nonanticipative investment strategies, in terms of the 
smallest positive solution to a parabolic partial differential inequality; 
this is determined entirely on the basis of the covariance structure of 
the model. The solution is intimately related to properties of strict 
local martingales and is used to generate the investment strategy 
which realizes the best possible arbitrage. Some extensions to non- 
Markovian situations are also presented. 

1. Introduction. In a Markovian model for an equity market with mean 
rates of return bi(3£(t)) and covariance rates ajj-(3£(i)), 1 < i,j < n, for its 
asset capitalizations X(t) = (X\(t), . . . ,X n (t))' G (0,oo) n at time t, what is 
the highest return on investment [as in (6.3) below] that can be achieved 
relative to the market on a given time-horizon [0, T], using nonanticipative 
investment strategies? What are the weights assigned to the different assets 
by such an investment strategy that accomplishes this? 

Answers: under suitable conditions, 1/U(T,X(0)) and 

Xi(t)Di log U(T - t, X(t)) + Xi{t) ^ {t) +Xn{t y i = 1, • • • , n, t e [0, T] , 

respectively. Here U : [0,oo) x (0,oo) n — > (0, 1] is the smallest nonnegative 
solution of the linear parabolic partial differential inequality 

dU 

(1.1) — (r,x)>£[/(r,x), (r,x)G(0,oo)x(0,oo) n , 
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subject to the initial condition £7(0, •) = 1, for the linear operator 

i=l j=l i=l \i=l n / 

with Di = d/dxi, D = (d/dxi, . . . ,d/dx n )' and Df- = d 2 /dxidxj. Further- 
more, U(T,X(0)) is the probability that the ([0,oo) n \{0})-valued diffusion 
process 2)(-) = (Yi(-), . . . , Y n (-))' with infinitesimal generator C as above and 
2)(0) = X(0) G (0,oo) n does not hit the boundary of the orthant [0,oo) n by 
time t = T. We note that the answers involve only the covariance structure 
of the market, not the actual rates of return; the only role these latter play 
is to ensure that the diffusion £(•) lives in (0,oo) n . 

Arbitrage relative to the market exists on [0,T], iff U(T,3Z(0)) < 1; this 
is deeply related to the importance of strict local martingales in the present 
context, and amounts to failure of uniqueness for the Cauchy problem 

dU 

— (r,x) =£{7(r,x), (r,x) G (0,oo) x (0,oo) n and U{0,-) = 1. 

Sufficient conditions for such failure of uniqueness are provided. 

Consider an "auxiliary market" with capitalizations 2)(-) = (Yi('),..., 
Y n (-))' as above. The probabilistic significance of the change of drift inherent 
in the definition of the operator £, from bj(x) for X(-) to X^=i(^.? a ij( x ))/ 
(xi + • • ■ + x n ) for 2)(-), is that it corresponds to a change of probability mea- 
sure which makes the weights fj(-) := Yi(-) / (Yi(-) -\ \- Y n (-)),i = 1, . . . ,n, 

of the auxiliary market portfolio martingales. Its financial significance is 
that it bestows to the auxiliary market portfolio i^(-) = (vi(-), ■ ■ ■ , ^n('))' the 
so-called numeraire property: any strategy's relative performance in the mar- 
ket with capitalizations 2)(-) is a supermartingale, so this market cannot be 
outperformed. This change need not come from a Girsanov-type (absolutely 
continuous) transformation; rather it corresponds to, and represents, the 
exit measure of Follmer (1972) for an appropriate supermartingale. 

Sections 2 and 3 set up the model, whereas Section 4 introduces the 
notion and offers examples of relative arbitrage; Section 5 makes the con- 
nection with strict local martingales. Section 6 formulates the problem, and 
Section 7 offers some preliminary results, actually in some modest generality 
(including non-Mar kovian cases). Section 8 sets up the Markovian model; 
the results are presented in earnest in Sections 9-11, Section 12 discusses a 
couple of examples in detail and a few open questions are raised in Section 
13. 

Related literature: the questions raised in this study are related to the 
work of Delbaen and Schachermayer (1995b). They bear an even closer con- 
nection with issues raised in the Finance literature under the general rubric 
of "bubbles" [see Definition 5 and Theorem 1 in Ruf (2009) for the precise 
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connection]. The literature on this topic is large, so let us mention the papers 
by Loewenstein and Willard (2000), Pal and Protter (2007) and, most sig- 
nificantly, Heston, Loewenstein and Willard (2007), as the closest in spirit to 
our approach here. We note the recent preprint by Hugonnier (2007), which 
demonstrates that arbitrage opportunities can arise in equilibrium models; 
this preprint, and Heston, Loewenstein and Willard (2007), can be consulted 
for an up-to-date survey of the literature on this subject and for some ex- 
plicit computations of trading strategies that lead to arbitrage. The need 
to consider state-price-density processes that are only local (as opposed to 
true) martingales has also been noticed in the context of "stochastic volatil- 
ity" models [e.g., Sin (1998), Wong and Heyde (2006)] and of pricing with 
long maturities [e.g., Hulley and Platen (2008)]. 

2. The model. We consider a model consisting of a money-market dB(t) = 
B(t)r(t) dt, B(0) = 1 and of n stocks with capitalizations, 

(2.1) dX i {t) = X i {t)(^ l (t)dt + ^a ik {t)dW k {t^ X l (0) = x i >0, 

for i = 1, . . . , n. These are defined on a probability space (f2, P) and are 
driven by the Brownian motion W(-) = (Wi(-), . . . ,Wk('))' whose K >n 
independent components are the model's "factors." 

We shall assume throughout that the interest rate process of the money- 
market is r(-) = 0, identically equal to zero; and that the vector-valued pro- 
cess £(•) = (Xi(-), . . . ,X n (-))' of capitalizations, the vector-valued process 
/?(•) = (/3i(-), . . . ,/3 n (-))' of mean rates of return for the various stocks and 
the (n x K) -matrix- valued process a(-) = {o~ik{-))i<i< n ,i<k<K of volatilities 
are all progressively measurable with respect to a right-continuous filtration 
F = {-7 7 (t)}o<t<oo which represents the "flow of information" in the market 
with ^"(0) = {0,17}, modP. Let a(-) := a(-)a' (■) be the covariance process 
of the stocks in the market, and impose for P-a.e. uj £ Vt the condition 

n „t 

(2.2) V/ (| / 8 i (t,w)|+Q ii (t,w))dt<oo VTg(0,oo). 
i=i J ° 

Under this condition the processes X\{-) , . . . , X n {-) can be expressed as 
Xi(-) = Xl exp{ / (/3i(t) - !<*«(*)) dt + Y%=x Jo <*ik(t) dW k (t)} > 0. 

In this setting, the Brownian motion W(-) need not be adapted to the 
"observations" filtration F. It is adapted, though, to the P- augmentation 
G = {G(t)}o<t<:oo of the filtration F, provided that K = n and that the 
matrix-valued process cr(-) is invertible — as in Assumption B below. 
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3. Strategies and portfolios. Consider now a small investor who decides, 
at each time t, which proportion 7r,(i) of current wealth V{t) to invest in the 
ith stock, i = 1, . . . , n; the proportion 1 — X^=i 7r «(*) =: n o(t) gets invested 
in the money market. Thus, the wealth V(-) = V v ' n (-) for an initial capital 
v G (0,oo) and an investment strategy ir(-) = (vri(-), . . . ,7r n (-))' satisfies the 
initial condition V(0) = v and 



dV{t) " M dXj(t) dB(t) 
(3.1) 

= ir'(t)[0(t)dt + (r(t) dW(t)\. 

We shall call investment strategy a G-progressively measurable process tt : [0, 
oo) x tt — > W 1 which satisfies for P-a.e. udJl the analogue 

cT 

(\ir'(t,L})/3(t,uj)\ +7C / (t,u})a(t,u))ir(t,u)))dt < oo, VTG (0,oo) 



o 

of (2.2). The collection of investment strategies will be denoted by %. 

A strategy ir(-) G % with Yh=i n(t,U}) = 1 for all (t,u) G [0,oo) x tt will 
be called portfolio. A portfolio never invests in the money market and never 
borrows from it. We shall say that a process 7r(-) is bounded, if for it there 
exists a real constant CV > such that ||7r(£,u;)|| < CV holds for all (t,uj) G 
[0,oo) x tt. We shall call long-only portfolio one that satisfies TTi(t,uj) > 
0, . . . ,7r n (t,ui) > 0,V(t,u;) G [0,oo) x tt, that is, never sells any stock short. 
Clearly, a long-only portfolio is also bounded. 

Corresponding to an investment strategy ir(-) and initial capital v > 0, 
the associated wealth process, that is, the solution of (3.1), is 

v v ' w (-) = wex p{^ ^(*) W) - ^r 7 ^)) dt + J Q 7r '(*) CJ (*) > o. 

The strategy g(-) = invests only in the money market at all times; it results 
in V v,e (-) = v, that is, in hoarding the initial wealth under the mattress. 

3.1. The market portfolio. An important long-only portfolio is the mar- 
ket portfolio; this invests in all stocks in proportion to their relative weights, 

(3.2) fr(t):=^Q, i = l,...,n,whereX(t):=X 1 (t) + ---+X n {t). 

Clearly V v ^{-) = vX(-)/X(0), and the resulting vector process fi(-) = (^i(-)> 
...,/i n (-))' of market weights takes values in the positive simplex A!L := 
{(mi,...,m n )' G (0,l) n | X£ =1 m* = 1} of R n . An application of Ito's rule 
gives, after some computation, the dynamics of this process as 



(3.3) d^ i {t) = n l {t) 



K 

^(t)dt + Y,rt k (t)dW k (t) 

k=l 



l,...,n. 
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Here r^(t) is the matrix with entries T p k (t) := <Tik(t) — Y^j=i MjW^jfcW) c « 
the ith unit vector in R n and the vector 7* i (t) := (7j t (t), . . . ,7n(f))' has 

(3.4) 7 f(*):=(e i -M(i)) / (/3W -«(%(*))• 

4. Relative arbitrage. The following notion was introduced in Fernholz 
(2002): given a real number T > and any two investment strategies ir(-) 
and /»(•), we call 7r(-) an arbitrage relative to p(-) over [0,T], if 

(4.1) P(V 1,ir (r) > V l ' p (T)) = 1 and F(V 1,7r (T) > V l,p (T)) > 0. 

We call such relative arbitrage strong if P(y 1 - 7r (T) > V" 1,P (T)) = 1. 

Arbitrage (resp., strong arbitrage) relative to q(-) = that invests only in 
the money market, is called just that, without the qualifier "relative." 



4.1. Examples of arbitrage relative to the market. Here are some exam- 
ples taken from the survey Fernholz and Karatzas (2009), especially Sections 
7 and 8, Remark 11.4, Examples 11.1 and 11.2. Suppose first that 



(4.2) ^ f x i (t)a ii (t)-^J2»i( t ) a v( t )»i( t )^ h V0<t<oo 

i=l i=l j=l 

holds almost surely for some constant h > 0. Then the long-only portfolio 
m(t) = m{t){c - log ia(t))/J(t),i = 1, . . . , n, J(t) := Y%=i Mj(*)(c - log fjtj(t)) 
is, for sufficiently large c> 0, a strong arbitrage relative to the market port- 
folio /i(-) over any time-horizon [0, T] with T> (21ogra)//i. 

Another condition guaranteeing the existence of strong arbitrage relative 
to the market is that there exists a real constant h > with 



(4.3) ^l(*)-"/*n(t) 



n - n n 

i=l i=l j=l 



>h 



V0<t<oo 



a.s. Then for c > sufficiently large, the long-only portfolio Tti(t) = X(t)(l/n) + 
(l-\(t))m(t), l<i<n, 1/A(i) :=l + (( f x 1 (t)---fji n (t)) 1 / n /c), is strong ar- 
bitrage relative to the market over any [0,T] with T> (2n 1 "( 1 / n ))//i. 



Remark 1. Suppose that all the eigenvalues of the covariance matrix- 
valued process a(-) are bounded away from both zero and infinity, uniformly 
on [0, oo) x $7, and that (4.2) holds. Then, for any given constant p £ (0, 1), 

the long-only portfolio /4 p) (t) = {f^i(t)) p (T,]=i(^j(t)) p )~ 1 , « = !,■•■,«, leads 
again to strong arbitrage relative to the market portfolio over sufficiently 
long time-horizons. It is also of great interest that appropriate modifications 
of the portfolio ^ p \-) yield such arbitrage over any time-horizon [0, T]. 
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5. Market price of risk and strict local martingales. We shall assume 
from now on that there exists a market price of risk $ : [0, oo) x Q — > M. K , an 
F-progressively measurable process that satisfies 

a(t,u))'8(t,cj) = /3(t,u>) V(t, uS) G [0, oo) X and 
(5-1) T 

Pf| \\$(t,oj)\\ 2 dt <oo, VTG (0,oo)J = 1. 

The existence of a market-price-of-risk process $(•) allows us to introduce 
an associated exponential local martingale, 

(5.2) Z(t) :=exp|-^ &'(s)dW(s) - ijf ||0(a) || 2 da J, 0<t<oo. 

This process is also a supermartingale; it is a martingale, if and only if 
K(Z(T)) = 1 holds for all T G (0,oo). For the purposes of this work it is 
important to allow such exponential processes to be strict local martingales; 
that is, not to exclude the possibility E(Z(T)) < 1 for some T E (0,oo). 
From (5.2) and (3.1), now written in the form 

(5.3) dV v ' n (t) = V v '^(ty(t)a(t)dW(t), W(t):=W{t)+ [ ${s)ds 

Jo 

on the strength of (5.1), the product rule of Ito's calculus shows that 

(5.4) Z (-)V v ' 7r (-) = v+ [ Z(t)V v ^(t)(a'(t)7T(t)-^t))'dW(t) 

Jo 

is a positive local martingale and a supermartingale, for every 7r(-) G T~L. 

If a(-) is invertible, we can take $(•) = cr'(-)a _1 (•)/?(•) as market price of 
risk in (5.1). If /?(•) = «(•)//(•) holds we can select #(•) = <r'(-)//(-) and get 
Z(-) = i'/y i; ' /i (-) = X(0)/X(-) from (5.4); there is then no arbitrage relative 
to the market because V v ' n {')/V v ' li (-) is a supermartingale for all tt(-) G 7^; 
thus E[V 1 ' 7r (T)/V 1 '' i (r)] < 1, a conclusion at odds with (4.1). 

5.1. Strict local martingales. Suppose the covariance process q(-) is 
bounded, and (4.1) holds for two bounded portfolios tt(-) and p(-). Then, 
for any market-price-of-risk process $(•) as in (5.1), the positive local mar- 
tingales Z(-) and Z(-)^' p (-) of (5.2), (5.4) are strict: E[Z(T)V V ' P (T)] < v, 
E(Z(T)) < 1 [Fernholz and Karatzas (2009), Section 6]. 

In particular, if the matrix a(-) is bounded, and (4.1) holds for some 
bounded portfolio ir(-) and for the market portfolio p(-) = /i(-) (these as- 
sumptions are satisfied, e.g., under the conditions in Remark 1), then 

E(Z(T))<1, E[Z(T)X(T)]<X(0), 

(5.5) 

E[Z(T)Xi(T)]<Xi(0), i = l,...,n. 
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6. Optimal arbitrage relative to the market. The possibility of strong 
arbitrage relative to the market, defined and exemplified in Section 4, raises 
an obvious question: what is the best possible arbitrage of this kind? 

One way to cast this question is as follows: on a given time-horizon [0,T], 
what is the smallest relative amount, 

(6.1) u(T):=mi{w>0\3ir(-)£Hs.t.V wX{0) ' n (T)>X(T), a.s.}, 

of initial capital, starting with which one can match or exceed at time t = T 
the market capitalization X(T)1 Clearly, < u(T) < 1; and for < w < 
u(T), no strategy starting with initial capital wX(0) can outperform the 
market almost surely, over the horizon [0, T]. That is, for every ir(-) G T~L and 
< w < u(T), we have p[y^(°)^(T) > X(T)] < 1. 

We shall impose from now on the following structural assumptions on the 
filtration F = {-7 r (t)}o<t<oo 5 the "flow of information" in the market. 

Assumption A. Every local martingale of the filtration F can be repre- 
sented as a stochastic integral, with respect to the driving Brownian motion 
W(-) in (2.1), of some G-progressively measurable integrand. 

Assumption B. We have K = n, and a(t) is invertible, Vi G [0,T]. 

Under these two assumptions, general results about hedging in so-called 
complete markets [e.g., Karatzas and Shreve (1998), Fernholz and Karatzas 
(2009), Section 10 or Ruf (2009)] based on martingale representation results, 
show that the quantity of (6.1) given as 

(6.2) u(T) = E[Z(T)X(T)]/X(p); that V U ^ X ^' 9 {T) = X(T) 
holds a.s. for some tt(-) G H; and that l/u(T) gives the highest return, 

(6.3) sup{g > 1| 3tt(-) G U s.t. V 1 ^ (T) > qV x ^{T), a.s.}, 

on investment, that one can achieve relative to the market over [0, T]. Arbi- 
trage relative to the market is possible on [0, T], if and only if u(T) < 1. 

The result in (6.2) provides no information about the strategy tt(-) that 
implements this "best possible" arbitrage, apart from ascertaining its ex- 
istence. In Section 8 we shall specialize the model of (2.1) to a Markovian 
context and describe tt(-) in terms of partial differential equations (Section 
11). We shall also characterize the quantity u(T) in terms of the smallest 
solution to a parabolic partial differential inequality, and as the probability 
of nonabsorption by time T for a suitable diffusion (Theorems 1, 2). 

Assumption A holds when F is (the augmentation of) ¥ w , the filtration 
generated by the Brownian motion VF(-); as well as when Assumption B 
holds, the /3j(-), cr^(-) are all progressively measurable with respect to = 
{.F*(£)}o<t<oo, ^(t) :=cj(X(s),0 < s<t), and F = F^ [Jacod (1977)]. 
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6.1. Generalized likelihood ratios. The positive local martingale Z(-)X(-), 
whose expectation appears in (6.2), can be expressed as 



* 1 



(6.4) Z(t)X(t)=X{0)-e X p{-J (m)'dW(s)-- J ||tf(s)|| 2 ds 
for < t < T. Here we have solved equation (5.4) for ir(-) = //(•) and set 

(6.5) £(•) := (•) - o-'(-)»(-), W(-) := W(-) + / 0(t) <ft, 

JO 

whence cr(-)^(') = /^(0 ~~ from (5.1); we thus re-cast (2.1) as 

(6.6) dXi(t)=Xi(t) 



V™ aij (t)Xj(t) " ~ 

x x (t) + ■ ■ ■ + x n {t) lkK ' feV ; 



On the other hand, we note from (6.4), (6.5) that the reciprocal of the 
exponential local martingale Z(-)X(-)/X(0) can be expressed as 

(6.7) A(-):=-^-=exp( f (&(t))' dW{t) -\ [ ||^)|| 2 dt); 



similarly, the reciprocal of the local martingale Z(-)Xj(-)/Xj(0) is 

(6-8) A,(.) := = ex P |jf'(^(t)) , dW?W(t) - 5 jf *}, 

where £«(.) :=-!?(■)- and W^(-) := W(-) + /o (*) 

Comparing (6.7) and (6.8), we observe that ^(O)A(-) = /ij(-)Aj(-) and cast 
the dynamics of (3.3) and (3.4) for the market portfolio //(•) as 



(6.9) d A t i (t) = A t i (t)(e i -At(t))V(t)dW(t) I i = l 



, . . . , 



n. 



If u(T) = 1, that is, Z(-)X(-) is a martingale on [0, T], no arbitrage relative 
to the market is possible on this time-horizon; the "reference" measure 

(6.10) F T (A):=E[Z(T)X(T)1 A ]/X(0), A G T(T), 

is a probability, that is, u(T) = Pt(^) = 1; and under Ft, the process 
W(t),0 <t < T, in (6.5) is a Brownian motion by the Girsanov theorem, 
so from (6.9) the market weights [J,i(t), . . . , fj, n (t),Q < t < T are martingales. 

We shall characterize next u(T) in terms of the Follmer exit measure, 
of a "generalized martingale measure" and of a measure Q with respect to 
which P is locally absolutely continuous [equations (7.3), (7.6)] and which 
plays, to a considerable extent, the role of Ft when Z{-)X{-) fails to be a 
P-martingale. The processes of (6.4)— (6.8) are important in this effort. 
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7. Exit measure of a positive supermartingale. We shall assume in this 
section that the process Z(-) of (5.2) is adapted to F = {.F(i)}o<t<oo and that 
this filtration is, in turn, the right-continuous version F{t) = f} £>0 J 70 (t + e) 
of a standard system ¥° = {J~°(t)}o<t<oo- to wit, each (fi, is isomor- 

phic to the Borel cr-algebra of some Polish space, and for any decreasing 
sequence {Aj}j & ^ such that Aj is an atom of F°(tj), for some increasing 
sequence {i/jjeN C [0, oo), we have fljeN A? ^ 0. 

The canonical example is the space f2 of right-continuous paths oj : [0, oo) — >• 
W 1 U {A}, where A is an additional "absorbing point" ; paths stay at A once 
they get there, that is, after T{oj) = inf{t > 0|cu(t) = A}, and are continu- 
ous on (0,T(w)). If F°{t) =a(u)(s),Q< s< t), then F° = {.F (i)}o<t<cx> is a 
standard system [see Follmer (1972), the Appendix]. 

Under these conditions, we can associate to the (P,F)-local martingale 
Z(-) • X(-) a positive measure ^ on the predictable cr-algebra of [0, oo] x Q, 

«p((T,oo]xA):=K[Z(T)X(T)l A ]/X(0), A £ F{T) , T £ [0, oo) , 

by invoking an extension result [Parthasarathy (1967), Theorem V.4.1, whence 
the assumptions on the nature of the probability space]. 

This is the "exit measure" of the supermartingale Z(-)X(-), introduced by 
Follmer (1972, 1973) [see also Delbaen and Schachermayer (1995a), Follmer 
and Gundel (2006)]. Follmer (1972) obtained a characterization of the (process- 
theoretic) properties of supermartingales, such as Z(-)X(-) here, in terms of 
the properties of ^3. It follows from his work that Z(-)X(-) is a: 

• martingale, if and only if in concentrated on {oo} x f2; 

• potential [i.e., u(oo) = 0], if and only if *p in concentrated on (0, oo) x f2. 

7.1. A representation of the Follmer measure. From Theorem 4 in Del- 
baen and Schachermayer (1995a) and Theorem 1 and Lemma 4 of Pal and 
Protter (2007), the process A(-) of (6.7) is a continuous martingale under 
some probability measure Q on the filtered space (£l,J-),F as above. The 
measure P is locally absolutely continuous with respect to Q, with ciP = 
A(T) ciQ on each J-(T); and the process W(-) of (6.5) is Q-Brownian motion 
[cf. Ruf (2009), Section 5]. Thus, from (6.9) the weights //i(-), ... ,//„(•) are 
martingales and satisfy X^iLiMiO) = 1 a - e -> under Q. 

We consider the first time the process A(-) hits the origin, 

(7.1) T := inf{t > 0|A(t) = 0} = inf{t > 0\Z(t)X(t) = oo} 

(infinite, if the set is empty). We have P(T < oo) = 0, but Q(T< oo) can 
be positive, so Q may not be absolutely continuous with respect to P; 
whereas, Q-a.e. on {T < oo}, we have Z(T + h)X(T + h) = oo, V/t > and 
Jo~ \ W(t)\\ 2 dt = oo. Intuitively, the role of the absorbing state A is to ac- 
count for events that have zero P-measure, but positive Q-measure. We also 
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introduce the first times the processes fJ-i(-) and Aj(-) hit the origin, 

(7.2) 77 := inf{t > 0|^(t) = 0}, % := inf{t > 0| Ai(t) = 0}. 



Proposition 1. (i) The quantity of (6.1) can be represented as 

(7.3) u(T)=*P((r,oo] xO) = Q(T>T). 

(ii) Suppose n>2 and that all capitalizations Xi(-), . . . ,X n (-) are real- 
valued Q-a.e. Then we also have the Q-a.e. representations 

(7.4) T = min Tf, as well as T = min % away from the event E, 

l<i<n l<i<n 

where E := {T < oo} n {^i(T) • • • fj, n (T) > 0}. TTtis event has Q-measure 
equal to zero, if for some real constant C > we have 

(7.5) \\^(t,u)\\ 2 <C{l + Tr(a(t,uj))) V(t,w) G [0, oo) x ft. 

Proof. We note <P((T,oo] x A) = E p (A- 1 (T)l An{r>T} ) = E Q (A(T) • 
A-^Tjl^n^rj) = Q(^n {T > T}), Wl G J"(T). With A = ft, we get (7.3). 
For A = {m{T) ■ ■ ■ n n (T) = 0}, this gives Q(^n{T > T}) = 0: all the weights 
[i\ (•))••• >/%(') are strictly positive [equivalently, all A"i (•),... ,X n (-) take 
values in (0,oo)] on [0,T), Q-a.e. 

Recall /ij(0)A(-) = /ij(-)Aj(-), Vi = l,...,n from (6.7), (6.8); this gives 
1 / A (-) = Er=i(^(°)/ A i(-)) on [0,T), and the first equation in (7.4). 

On the event {T < oo} \ £7, for some j G {1, . . . , n} we shall have ^j(T) = 
0, thus also Tj = T and the second equation in (7.4). On the other hand, 
we have seen that % = oo, Vi = 1, . . . ,n holds Q-a.e. on {T = oo}, so this 
equation is valid on {T = oo}. 

Finally, from (6.6), (6.7): Tr(a(t,w)) eft < oo, J T ||^(t,o;)|| 2 (it = oo for 
Q-a.e. w G £ C {T < oo}. Then (7.5) implies J T ||i?(t, w)|| 2 dt < oo, and 
= #(•) - o-'(0m(-) gives J r ||tf(t,a/)|| 2 cft<oo, thus Q(£) = 0. □ 

Equation (7.3) can be thought of as a "generalized Wald identity" [cf. 
Problem 3.5.7 in Karatzas and Shreve (1991)]. In Section 9.3 we shall obtain 
a characterization of the type (7.3) in a Markovian context, in terms of 
properties of an auxiliary diffusion and with the help of an appropriate 
partial differential equation. This will enable us to describe the investment 
strategy that realizes the optimal arbitrage. 
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7.2. A generalized martingale measure. In a similar vein, there exists on 
the filtered space (0, J-),¥ a probability measure Q under which 

L(t):=l/Z(t) = exp(J^(s)dW(a)-^J \\${s)\\ 2 da J, 0<t<oo, 

is a martingale, and dF = L(T) dQ on each J~(T), whereas W(-) of (5.3) is 
Q-Brownian motion. Under Q, the processes Xi(-), i = 1, . . . , n are nonnega- 
tive local (and super-) martingales, dXi(t) = Xi(t) J2k=i °~ik{t) dW k (t) . This 
justifies the appellation "generalized martingale measure" for Q. 

Defining S := inf{i > 0\L(t) = 0}, we have P(«S < oo) = and Z(-) is a 
strict P-local martingale if and only if Q(<S < oo) > [a potential, if and 
only if <Q)(<S < oo) = 1]; and the expression of (6.1), (6.2) is 

(7.6) u(T)=e\x(T)/X(0))1 {s>t} ]. 

This last expression takes the formu(T) = 1-E%X(S)/X(0))1 {(S < T} ] when 
X(-AT) is a Q-martingale; from (5.3), this will be the case under the Novikov 

condition E®[exp{i fi'(t)a(t)n(t) dt}] < oo. Moreover, u(T) = 1 (no arbi- 
trage relative to the market is possible on [0,T]), if and only if: X(- A T) is 
a Q-martingale, and X(S)1{s<t} = holds Q-a.e. 

8. A diffusion model. We shall assume from now on that K = n and that 
the processes /%(•), o~i k (-), 1 < i, k < n in (2.1) are of the form 

(8.1) pi(t)=bi(X(t)), a ik (t)=s ik (X(t)), 0<t<oo. 

Here X(t) = (X\ (t), . . . , X n (t))' is the vector of capitalizations at time t, and 
bj : (0, oo) n — > E, Sik ■ (0, oo) n — > M are continuous functions. We shall denote 
by b(-) = (bi(-), • ■ - ,*>„(•))' and s(-) = (s ik (-))i<i<n,i<k<n the vector and ma- 
trix, respectively, of these local rate-of-return and local volatility functions. 
With this setup, the vector process 3t(i),0 < t < oo of capitalizations be- 
comes a diffusion, with values in (0,00)™ and dynamics 

n 

(8.2) dX l (t) = b l {3i{t))dt + Y,Sik(X(t))dW k (t), i = l,...,n, 

k=l 

where for x (xi,...,x n )' G (0,oo) n we set a^(x) := YZ=i s *fc(x)s J) - fc (x), 

(8.3) bj(x) := Xibj(x), S ifc (x) := XjS ifc (x), Oij(x) := a^a^x). 
This diffusion X(-) has infinitesimal generator 

1 n n 71 

(8-4) £/ := - E M*)^/ + E b i( x ) D if- 

i=l j=l i=l 
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Assumption C. For every x e (0,oo) n , the matrix s(x) = (s i: ,(x))i<ij< n 
is invertible; the system (8.2) has a unique-in-distribution weak solution, 
with 3t(0) = x and values in (0,oo) n ; and for 0(x) := s _1 (x)b(x), the fol- 
lowing analogue of (2.2), (5.1) holds for each T E (0,oo): 

n ,. T 

(8.5) V/ (MX(t))\+aii(X(t)) + e](X{t)))dt<oo a.s. 

It follows from this assumption that the Brownian motion W(-) is adapted 
to the augmentation of the filtration F , and that ■$(•) = &(£(■)) is a market- 
price of risk process as postulated in (5.1). The following conditions from 
Bass and Perkins (2003), in particular their Theorem 1.2 and Corollary 1.3, 
are sufficient for the existence of a weak solution for (8.2) which is unique in 
distribution: the functions Sjfc(-), bj(-) of (8.3) can be extended by continuity 
on all of [0,oo) n ; bj(-) and f)ij(x) := y^a^TTa^- (x) > are Holder continuous 
on compact subsets of [0,oo) n ; and we have 

bj(x) > for Xi = 0; 

(8.6) 

||b(x)|| + ||s(x)|| < C(l+ ||x||) VxG [0,oo) n , 
and f)ij(x) = for j, x S O n , where O n is the boundary of [0, oo) n . 

Remark 2. The diffusion £(•) of (8.2) takes values in (0,oo) n , if and 
only if the diffusion S(-) = (Hi(-), . . . , S n (-))', Hf(-) := 1/Xj(-), with dynamics 

n 

(8.7) dHi(t) = q i (H(t))dt + 5^t ifc (H(t))dWfc(t) J i = l,...,n, 

k=l 

and t ifc (£) := -^s ifc (l/Ci, • • • , 1/Cn), q%(Q ■= Ci(^i-bi)(l/Cl, • • • , takes 
values in (0,oo) n . Thus, any conditions guaranteeing the existence of a non- 
explosive solution to the SDEs of (8.7) for all times, such as linear growth 
for cjj(-) and tjfc(-), also ensure that X(-) takes values in (0,oo) n . 

Alternatively, one may invoke results of Friedman (2006), Section 9.4 and 
Chapter 11, to obtain conditions on bj(-),Sifc(-) under which the diffusion 
£(•) of (8.2) never attains any of the faces {x\ = 0}, . . . , {x n = 0} of O n . In 
particular, if these functions can be extended by continuity on all of [0,oo) n ; 
the Sjfc(-) are continuously differentiable; the matrix o(-) degenerates on the 
faces of the orthant; and the so-called Fichera drifts 

1 n 

(8.8) f i (x):=b i (x)--£z> i a ii (x) 

3=1 

are nonnegative on {x^ = 0}, for each i = 1, . . . ,n; then £(•) takes values in 
(0,oo) n [see Friedman (2006), Theorem 9.4.1 and Corollary 9.4.2]. 
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Assumption D. There exists H : (0, oo) n — > R of class C 2 , such that 

(8.9) b(x) = a(x)DH(x), Vxe(0,oo) n . 

In light of Assumption C, this new requirement amounts essentially to 
postulating that the vector field a -1 (-)b(-) be conservative; it is imposed here 
for technical reasons (cf. discussion in Remark 3). Under it, the generator of 
(8.4) becomes £/(x) = £? =1 £"=i a ij (x)[i£>?./(x) + A/(x)D^(x)], and 
we have 

(8.10) 6(x)=s'(x)M(x) and s(x)6(x) = b(x), xG(0,oo) n . 

Throughout the remainder, Assumptions B, C, D will be in force, and 
F = Ff ; this is a natural choice, and consistent with Assumption A. 

9. A parabolic PDE for the function U(t, x). The uniqueness in dis- 
tribution posited in Assumption C implies that X(-) is strongly Marko- 
vian; we shall denote by P x the distribution of this process started at 
X(0) = x G (0,oo) n . Our objective now is to study 

(9.1) U(T,x) := E F *[Z(T)X(T)]/( Xl + ■■■ + x n ), 

the quantity of (6.1), (6.2) in this diffusion context. We start by observing 
that with -ff(-) as in Assumption D and the notation of (8.4) and (8.10), 
Ito's rule gives H(X(T)) - H(X(0)) - f*£H(X(t))dt = &' '(£(*)) dW \t) , 
and the exponential local martingale Z{-) of (5.2) becomes 

(9.2) Z(-) = expjtf(£(0)) - H (%(•)) -J k(X(t))dt\. 
In particular, Z(-) is F-^-adapted. We are setting here 



n n 



(9.3) fc(x) := - ^PJ^W + A^(x)^if(x)], 

i=i j=i 

n 

(9.4) g(x) := e ^( x )^xi, G(T,x) := E FX [g{3L{T))e~ ^ k ( x( - t))dt }. 

i=l 

With this notation, the function of (9.1) becomes U(T,x) = G(T,x)/g(x). 
A bit more generally, these considerations — coupled with the Markov prop- 
erty of £(•) — lead for any < t < T to the a.s. identity 



¥F~ [X {T)Z{T)\F{t)] _G(T-t,y) 



X{t)Z{t) g(y) 



= U{T-t,X(t)). 



The following Assumption E will also be imposed from now onward. It 
amounts to assuming that the function U(-, •) of (9.1) is of class C 1,2 . Note 
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that (9.6) is satisfied, at least in the support of £(•), thanks to the assump- 
tion [/(•, •) E C 1,2 ((0,oo) x (0,oo) n ) and to the P x -martingale property of the 
process G{T - t,X(t))e~ti fc (*M) d «,0 < t < T. 

Assumption E. The function G(-,-) in (9.4) takes values in (0, oo), is 
continuous on [0, oo) x (0,oo) n , of class C 1 ' 2 on (0, oo) x (0, oo) n , and solves 

BC 

(9.6) — (t,x)=£G(t,x)-£:(x)G(t,x), r € (0,oo),x € (0,oo) n , 

(9.7) G(0,x)= 5 (x), xG(0,oo) n . 

This Cauchy problem is exactly the one arising in classical Feynman-Kac 
theory [see, for instance, Friedman (2006), Sections 5.6, 6.5, Karatzas and 
Shreve (1991), Section 5.7 and Janson and Tysk (2006)]. From Theorem 
1 and the remark following it in Heath and Schweizer (2000), Assumption 
E holds if: the functions bj(-), Sjfc(-) of (8.3) are continuously differentiable 
on (0,oo) and satisfy the growth condition in (8.6); the functions %•(■) of 

(8.3) satisfy the nondegeneracy condition (9.14) below; the function g(-) in 

(9.4) is Holder continuous, uniformly on compact subsets of (0,oo) n ; the 
continuous function k(-) of (9.3) is bounded from below; and the function 
£?(•,•) in (9.4) is real-valued and continuous on (0, oo) x (0,oo) n . This latter 
requirement is satisfied, for instance, if the functions tifc(-)> C|j(-) in (8.7) 
obey linear growth conditions, and the function g(£) := g(l/£i, ■ ■ ■ , l/£n) 
has polynomial growth [see Karatzas and Shreve (1991), Problem 5.3.15, as 
well as Heath and Schweizer (2000), Lemma 2 (and the paragraph preceding 
it)]. 

Sustained computation shows then that the Cauchy problem of (9.6), (9.7) 
for G (•,•), leads to a corresponding Cauchy problem for {/(-,•), namely 

dU , , 1 AA , . . , r^v^ fliiM-DjE/Yr, x) 

<»■»> *) = 2 £ £ °«<*>^, x >+EE £ +. 

j=l 3=1 i=l 3=1 

for (r,x) e (0,oo) x (0,oo) n ; and 17(0, x) = 1 for x £ (0, oo) n . 

9.1. An informal derivation of (9.8). Rather than including the compu- 
tations which lead from (9.6) to equation (9.8), we present here a rather 
simple, informal argument that we shall find useful also in the next subsec- 
tion, in a more formal setting. We start by casting (6.4) as 

d{ x(^z(t) = E ( E * (t)CTjfc ( * } " ^ fc (t) ) dw/fc (t) = - E ®* (*(*)) dW * (*) > 
fe=i \ t=i / fe=i 
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where, by analogy with (6.5), we have set 

(9.9) e fcW :=6 iW -t(^L), 

1 = 1 
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k = 1, . . . , n. 



On the other hand, assuming that U(-, •) of (9.1) is of class C > , we obtain 
from Ito's rule and with R k (r, x) := J2?=i x i s ik( x )DiU(T,x), k = 1, . . . , n, 

(Pitt \ n 
C-U — — — J (T — i, £(*)) dt + ^ 2? fc (T - i, £(*)) dWfe(t). 
r ' fe=i 

The product rule of the stochastic calculus applied to the process 
(9.10) N(t) := X(t)Z(t)U(T - t, X(t)) = ^{X{T)Z{T)\F{t)} 
of (9.5), leads then to 



dN(t) 
X(t)Z(t) 



dU(T - t, X(t)) + U(T - t, X(t)) 



d(X(t)Z(t)) 



X(t)Z{t) 

n 
k=l 

C(T-t,X{t))dt 

n 

+ ^[i? fc (r - 1, x(t)) - u(t - 1, x(t))e k (x(t))} dw k (t). 



k=i 



We have set 



,x):=(cU- ?pj (t, Mr, x)0 fc (x) 

^ T ' k=i 



- £ £ ay (x)Z^ Z7(r, x) + £ £ ^ ^ . - ^ (r, x) , 



i=l j=l 



i=l 3=1 



xi H Yx r 



where the last equality is checked easily from (8.4) and (8.10). But the 
process N(-) of (9.10) is a martingale, so the term C(r,x) should vanish, 
and 



(9.11) 



dN(t) 



■E 
fc=i 



R k (T-t,X(t)) 
U(T-t,X(t)) 



@ k (X(t)) 



dW k (t). 



In other words, the function E/ : [0, oo) x (0, oo) n — > (0, 1] of (9.3) must satisfy 
the parabolic partial differential equation (9.8), as postulated earlier. 
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Remark 3 . This informal derivation suggests that it may be possible to 
dispense with Assumptions D, E altogether, if it can be shown from first prin- 
ciples that the function U of (9.1) is of class C 1 ' 2 on (0, oo) x (0,oo) n . Indeed, 
under suitable conditions, one can rely on techniques from the Malliavin 
calculus and the Hormander hypoellipticity theorem [Nualart (1995), pages 
99-124] to show that the (re + 2)-dimensional vector (X(T),T(T),E(T)) with 

T(T) := J^Q(X(t)YdW(t) and 3(T) := J Q T \\Q(X(t))\\ 2 dt has an infinitely 
differentiable probability density function, for any given T G (0,oo). This 
provides the requisite smoothness for the function 

E/( T ,x) = E px [(X!(T) + • • • + X n (T))e r(T ^ E ^^}. 

x\ H h x n 

The conditions needed for this approach to work are strong; they include 
the infinite differentiability of the functions Sik(-), 1 < i,k < n, as 

well as additional algebraic conditions which, in the present context, are 
somewhat opaque and not very easy to state or verify. For these reasons we 
have opted for sticking with Assumptions D, E; these are satisfied in the 
Examples of Section 12, are easy to test and allow us to represent Follmer's 
exit measure via (9.23), (9.24) without involving stochastic integrals. 

9.2. Results and ramifications. Equation (9.8) is determined entirely from 
the volatility structure of model (2.1). Furthermore, the Cauchy problem of 
(9.8), U(0, •) = 1, admits the trivial solution f7(r,x) = 1; thus, the existence 
of arbitrage relative to the market portfolio over a finite time-horizon [0,T] 
is tantamount to failure of uniqueness for the Cauchy problem of (9.8), 
U(0, •) = 1 over the strip [0,T] x (0,oo) n . 

Remark 4. Assume there exists some h > such that the continuous 
functions ajj(-), 1 < i,j < n satisfy either of the conditions 

n n n 

(9.12) (xH hx n )y^^a„(x) - (x) > h(x 1 H h x n ) , 



(9.13) {x x - 



n ) 



i i=i j=i 

n 1 n n 

1 !>(*)- 

.i=l i=l j=l 



> h(xi H h x n ) 



for all x S (0, oo) n [we have just re-written (4.2) and (4.3) in the present 
context]. Then from the results reviewed in Section 4 we deduce that, for all 
T > (21ogn)//i under (9.12), and for all T > (2n 1 ~( 1 / n ))//i under (9.13), we 
have C/(T,x) < 1, VxG (0,oo) n . In particular, under either (9.12) or (9.13), 
uniqueness fails for the Cauchy problem of (9.8), U (0, •) = 1. 

Whenever uniqueness fails for this problem, it is important to know how 
to pick the "right" solution from among all possible solutions, the one which 
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gives the quantity of (9.1). The next result addresses this issue; it implies 
that G(-, •) in (9.4) is the smallest nonnegative, continuous function, of class 
C 1 ' 2 ((0,oo) x (0,oo) n ), which satisfies (dG/dr) > CG - kG and (9.7) [cf. 
Karatzas and Shreve (1991), Exercise 4.4.7 for a similar situation]. 

Theorem 1. The function U : [0, oo) x (0,oo) n — > (0,1] of (9.1) is the 
smallest nonnegative continuous function, of class C 1 ' 2 on (0,oo) x (0,oo) n , 
that satisfies U(0,-) = 1 and (1.1). 

Proof. Consider any continuous function U : [0,oo) x (0,oo) n — > [0,oo) 
which is of class C 1 ' 2 on (0, oo) x (0,oo) n , and satisfies (1.1) and U(0, ■) = 1 
on (0,oo) n ; we shall denote by il the collection of all such functions. We 
introduce N(t) := X(t)Z(t)U(T - t,3L(t)), < t < T as in (9.10). 

Repeating verbatim the arguments in Section 9.1, we use (1.1) to conclude 
that the nonnegative process N{-) is a local supermartingale. Thus N(-) is 

bona-fide supermartingale, (a* + h x n )U(T,x) = iV(0) > E px (iV(T)) = 

E px (X(T)Z(T)) holds for every (T, x) £ (0,oo) x (0,oo) n , and U(T,x) > 
U(T,x) follows from (9.1). □ 

Proposition 2. Assume that the continuous functions (cijj(-))i<ij<n of 
(8.3) satisfy the following nondegeneracy condition: for every compact subset 
fC of (0, oo) n ; there exists a number e = ex. > such that 

n n 

(9.14) ( z )&& ^ e ll£H 2 > Vz € /C,e £ M n - 
i=i j=i 

Then, if 

(9.15) U(T, x) < 1 for some x £ (0, oo) n 
holds for some T £ (0,oo), we have 

(9.16) U(T, x) < 1, V(T, x) £ (0, oo) x (0, oo) n . 

Proof. Let us work first under the stronger assumption 

(9.17) f/(T,x)<l, Vx£(0,oo) n , 

for some T £ (0,oo). For every r > 0, we consider the set <S(r) := {x £ 
(0,oo) n |C/(r,x) = 1} and define r* := sup{r £ (0,oo)|5(r) / 0} (with r* = 
if the set is empty). Assumption (9.17) amounts to r* < oo, and the claim 
(9.16) to r* = 0; we shall prove this claim by contradiction. 

Suppose > 0; then C/(r* — 5, x*) = 1 for any given 5 £ (0,t*/2), and 
some x* £ (0,oo) n . For any given x £ (0,oo) n , consider an open, connected 
set D which contains both x and x*, and whose closure D is a compact 
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subset of (0,oo) n ; in particular, we have inf { ||y — z|||z G D, y G O n } > 0. 
The function U(-,-) attains its maximum value over the cylindrical do- 
main G; = {(r,£)|0 < r < t* + 1,£ G .D} at the point (r* — 5,x*), which lies 
in the interior of this domain. By assumption then, the operator Cf = 
(1/2) ELi £j=i + E"=i^(x)A/ of (1.2) with 

(9.18) £;(*):= Wb^x), b t (x):=f; , i = l,...,„, 

is uniformly parabolic with bounded, continuous coefficients on l£, so from 
the maximum principle for parabolic operators [Friedman (2006), Chapter 
6], 

(9.19) U(r, x) = 1 V(r, x) G [0, t* - 5) x (0, oo) n . 

Now let us recall the P x -a.s. equality E px [X(T)Z(T)|J r (t)] = U(T-t,X(t)) • 
from (9.5); we apply it with < t < r* - 5, < T - t < r* - 5, then 
take expectations with respect to the probability measure P x , and use (9.19) 
along with (9.1), to obtain for every T G [0, 2(t* — 5)], 

But since 2(t* — 5) > r*, this contradicts the definition of 

Now we revert to (9.15); as J. Ruf (private communication) observes, yet 
another application of the maximum principle, as above, leads to (9.17). □ 

Corollary. Under the nondegeneracy condition (9.14), and with either 
(9.12) or (9.13), inequality (9.16) holds. That is, arbitrage with respect to 
the market exists then over any time-horizon [0,T] with TG (0, oo). 

9.3. An auxiliary diffusion. Let us consider now the diffusion process 
2)(-) with infinitesimal generator C as in (1.2), (9.18) and dynamics 

n 

(9.20) dY i (t)=%(%(t))dt + ^s ik (m))dW k (t), i = l,...,n. 

k=\ 

Assumption F. The system of SDEs (9.20) admits a unique-in-distribution 
weak solution with values in [0,oo) n \{0}. 

This will be the case, for instance, if the drift functions bj(-), 1 < i < n 
of (9.18) can be extended by continuity on all of [0,oo) n and satisfy the 
Bass and Perkins (2003) conditions preceding, following and including (8.6). 
The resulting process 2)(-) is then Markovian, and we shall denote by Q y 
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its distribution with 2)(0) = y 6 [0,oo) n . Unlike the original process £(•), 
which takes values in (0,oo) n , this new process 2)(-) is only guaranteed to 
take values in the nonnegative orthant [0,oo) n \ {0}. In particular, with 
x G (0, oo) n the first hitting time 



of the boundary O n of [0, oo) n may be finite with positive Q x -probability. 

Our next result shows that this possibility amounts to the existence of 
arbitrage relative to the market, and to the lack of uniqueness for the Cauchy 
problem of (9.8) and [7(0, •) = 1. 

Theorem 2. With the above notation and assumptions, including (9.14), 
the function U : [0, oo) x (0,oo) n — > (0, 1] of (9.1) admits the representation 



Proof. The function on the right-hand side of (9.22) is space-time har- 
monic for the diffusion 2)(-) on (0,oo) x (0,oo) n , so it solves equation (9.8) 
there [cf. Janson and Tysk (2006), Theorem 2.7]. Consider any function V 
in the collection il of Theorem 1; then V(T — t, 2)(t))l{£>t}> < t < T is a 
nonnegative local (thus a true) Q x -supermartingale, and we deduce 

V(T, X ) >E« x [F(0,2)(T))l {x>r} ] = Q X (T>T), 

(T,x) G (0,oo) x (0,oo) n . 

The claim follows now from the proof of Theorem 1. □ 

Corollary. Under the assumptions of Theorem 2, for any given x G 
(0,oo) n the W~ -supermartingale Z(-)X(-) is under a: 

• martingale, if and only if Q X (T < oo) = 0; 

• potential [i.e., lim T _, 0O E px (Z(T)X(r)) = 0/, Q x (1 < oo) = 1; 

• strict local (and super-) martingale on any time-horizon [0,T] with TG 
(0,oo), if and only if Q X (T < oo) > 0. 

We represent by analogy with (7.3) the exit measure *}3 X of the super- 
martingale Z(-)X(-) with initial configuration 3t(0) =x, in the form 



(9.21) 



1:=mi{t >0|2)(t) GO n } 



(9.22) 



f/(T,x) =Q X [T>T] 



(T,x) G (0,oo) x (0,oo) n . 



(9.23) <p x ((T,oo] x n) =U(T,x) =Q*[T>T] 



and from (9.2)-(9.5) we have for A G T(t),0< t < T, 




e -J*k(X(8))d8 
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When x £ (0,oo) n and the quantity of (9.22) is equal to one, the Q x - 
distribution of the process 2)(i),0 < t < T in (9.20) is the same as the P x - 
distribution of the original stock-price process X(t),0 <t<T; this follows by 
comparing (9.20) and (9.18) with (6.6), and denoting by Pj. the probability 
measure Pt of (6.10) with X(0) = x. We have in this spirit the following 
result, by analogy with Remark 2. 

Proposition 3. Under the assumptions of Theorem 2, suppose that the 
functions Sjfc(-) are continuously differentiable on (0,oo) n ; that the matrix 
a(-) degenerates on O n ; and that the analogues of (8.8), the Fichera drifts 

(9.25) f,(x) := b*(x) - | =E( - " tj(X | T " 5 W*)) 

z J=1 j=iV Xl_l l_Xri / 

/or the process 2)(-) o/ (9.20), can be extended by continuity on [0, oo) n . If 
fi( - ) > ZioZds on eac/i /ace {xj = 0}, i = 1, . . . , n of the orthant, then we have 
[/(•,•) = 1 in (9.22), and no arbitrage with respect to the market portfolio 
exists on any time-horizon. 

If, on the other hand, we have fj(-) < on each face {xi = 0} of the 
orthant, then U(-,-) < 1 in (9.22) and arbitrage with respect to the market 
portfolio exists, on every time-horizon [0,T] with T £ (0,oo). 

Proof. In light of Theorem 2, the first claim follows from Theorem 
9.4.1, Corollary 9.4.2 of Friedman (2006), and the second is a consequence 
of the support theorem for diffusions [Ikeda and Watanabe (1989), Section 
VI.8]. □ 

Remark 5. (i) The "relative weights" ^(t) := Y i (t)/(Y 1 (t) H \-Y n (t)), 

i = 1, ... ,n have dynamics similar to (6.9), 

(9.26) dui{t) = ui(t)(ti - i/(t))'s(g)(t)) dW{t). 

They are thus Q x -martingales with values in [0, 1] (cf. Section 6.1); so, when 
any one of them hits either boundary point of the unit interval, it gets 
absorbed there. In terms of them, the first hitting time of (9.21) can be 
expressed as in (7.4), T = mhii<j< n Tj, where Tj := inf{t > 0|z/j(£) = 0}. 

(ii) The measure Q x corresponds to a change of drift, from b(-) in (8.2) 
to b(-) in (9.18), (9.20); this ensures that, under Q x , the components of 
the new, "fictitious" market portfolio are martingales, that v(-) has the 
numeraire property, and thus that v{-) cannot be outperformed. 
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10. Markovian market weights. Let us assume now the form 

bj(x) = *8 i (x 1 /x,...,x n /x), Sjfc(x) = & ik (x 1 /x,...,x n /x) 

for the functions of (8.1), with x := Y^=i x i an d suitable continuous func- 
tions 2$i(-)> o n A™. For m= (mi,...,m„)' G A™, we set ^4jj(m) := 
E&=i ®ifc( m ) ( 3jfc( m )- In words, we consider instantaneous growth rates and 
volatilities that depend at time t only on the current configuration fj,(t) = 
(fj,i(t), . . . , Hn(t))' of relative market weights, so the process //(•) of (3.3) is 
now a diffusion with values in the positive simplex A" and 



(io.i) din(t) = m(t) 



T l (p(t))dt + Y / TiMt))dW k (t) 



fc=l 



l,...,ra, 



with Tik(m) := 6 ifc (m) - £ " =1 m,e^(m),^(m) := ££ = l^(m)7^(m), 

n n n n 

Tj(m) := 5Sj(m) — ^~^mj*Bj(m) — ^~]mjAjj(ia) + ^^mjm^(m). 

3=1 3=1 3=1 fe=l 

In this setup, the function of (9.1) can be expressed in the form £7(T, x) = 
Q(T, xi/x, . . . , x n /x), in terms of a function Q:(0,oo) x A" — > (0, 1] that 
satisfies the initial condition Q(0, •) = 1 and the equation 

— (r,m) = - ^2^2m l m j V ij (m)Df j Q(r,m), (r,m) G (0,oo) x A£, 
r i=i 3 -=i 

which appears on page 56 of Fernholz (2002) and can be derived from (9.8). 
On the other hand, by analogy with Theorem 2 and (9.26), the quantity 
Q(T,m) is the probability that the process u{-) = (y\ (•),..., v n ('))' with 
z^(0) = m G A" and dynamics (10.2) below, does not hit the boundary of 
the nonnegative simplex A n := {m G [0, l] n | Ei=i m i = x l before t = T: 

n 

(10.2) dv.it) = v l {t)Y,Tik{v{t))dW k {t), i = l,...,n. 

k=l 

11. The investment strategy. Let us substitute now the expressions of 
(9.9) into (9.11), to obtain the dynamics of the martingale N(-) = Z(-)X{-)U(T- 
■,£(•)) in (9.10), with JV(0) =£:=X(0)U(T,X(0)), 

n rt 

N (f) = Z + J2 N(s)^ k (T-s,X(s))dW k (s), 0<t<T, 
k=i Jo 

* fc (r,x) := Vs ifc (x)fxiAlog?7(r,x) + ) - 6 fc (x). 
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Thus we can identify the "replicating strategy" tt(-) of (6.2) as 

(11.1) 9 i (t) = X i (t)D i ]ogU(T-t,X(t)) + (X i (t)/X(t)), i = l,...,n, 

and its value as V^(t) = N(t)/Z(t) = X(t)U(T - t, £(*)), < t < T. 

Remark 6. In the special case of a Markovian model (10.1) for the 
market weights of /i(-) = (fJ>i(-), ■ ■ ■ ,fJ, n (-))', expression (11.1) takes the form 

7Vi(t) = fn(t) ^1 + A log Q(T - t, fi(t)) - J2 H(t) D 3 l °S Q(T - t, fi(t))j 

of a "functionally-generated portfolio" in the terminology of Fernholz (2002), 
page 56; whereas the value is V*' 9 (t) = X(t)Q(T - t,fi(t)), < t < T. 

In this case we have Y27=i ^i(') = 1 ; the strategy that implements the best 
possible arbitrage relative to the equity market never borrows or lends. 

12. Examples. We discuss in this section two illustrative examples. Addi- 
tional examples, in which the investment strategy tt(-) of (11.1) that realizes 
the optimal arbitrage can be computed in closed form in dimension n = 1, 
can be found in Ruf (2009). 

For the first of these examples, take n = 1, [3(t) = 1/X 2 (t) and a(t) = 
l/X(t) in (2.1) where the process X(-) satisfies dX(t) = (1 / X (t)) dt + dW (t) 
and X(0) = 1. This is a Bessel process in dimension three — the radial part of 
a 3-D Brownian motion started at unit distance from the origin — and takes 
values in (0,oo). We have then #(t) = 1/X(t), Z{t) = l/X{t) for < t < oo in 
(5.1) and (5.2), so Z(-)X(-) is very clearly a martingale. However, Z(-) is the 
prototypical example of a strict local martingale — we have E(Z(T)) < 1 for 
every T G (0,oo) [e.g., Karatzas and Shreve (1991), Exercise 3.36, page 168]. 
This example is taken from Karatzas and Kardaras [(2007), page 469], where 
an arbitrage with respect to the money-market is constructed in closed form. 
It illustrates that it is possible for Z{-) to be a strict local martingale and 
Z{-)X{-) to be a martingale; in other words, the second and third inequalities 
in (5.5) fail, while the first stands. 

Here we have G(x) = 1/x, H(x) = logx and k(-) = 0, g(-) = 1, G(-, •) = 1 in 
(9.3), (9.4), thus U(T, x) = 1 for all T G [0, oo), x G (0, oo). Arbitrage relative 
to X(-) does not exist here, despite the existence of arbitrage relative to the 
money market and the fact that Z(-) is a strict local martingale. Note that 
b(x) = 1/x in (9.18), so the diffusion of (9.20) is again a Bessel process in 
dimension three, dY(t) = (1/Y(t))dt + dW(t), Y(0) = y > 0. This process 
never hits the origin, so the probability in (9.22) is equal to one, for all 
TG [0,oo). 
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12.1. The volatility-stabilized model. Our second example is the model 
of "stabilization by volatility" introduced in Fernholz and Karatzas (2005) 
and studied further by Goia (2009). With n > 2, £ G [0, 1] this posits 

&(*) = (l + C)/(2/Zi(t)), 

(12.1) 

o-jfe(t) = 6 ik (m(t)) 1/2 ; l<i,fc<n, 

that is, rates of return and volatilities which are large for the small stocks and 
small for the large stocks. The conditions of Bass and Perkins (2003) hold 
for the resulting system of SDEs in the notation of (3.2) with k := (1 + Q/2, 



(12.2) dXi(t) = nX(t) dt + y/Xi(t)X(t) dWi(t), i = l, 



,n. 



The unique-in-distribution solution of (12.2) is expressed in terms of in- 
dependent Bessel processes JHi(-), . . . ,9t n (-) in dimension An with Xi(t) = 
mf{A(t)) > and A(t) := (1/4) f*X(s)ds. In particular, £(•) takes values 
in (0,00)"; for more details on these Lamperti-like descriptions and their 
implications, see Fernholz and Karatzas (2005) and Goia (2009). Condition 
(8.5) is satisfied in this example, so Assumption C also holds. 

For the model of (12.1), we have Oj(x)/k = su(x) = ((#iH Vx n ) jxi) 1 ! 2 , 

b i (x) = K(x 1 -\ \-X n ), t)ij(x) = S i j(x 1 -\ \-X n ), Oy(x) = Xit)ij(x) 

for 1 < i,j < n. The assumptions of Theorem 2 and of Propositions 1 and 
3 are all satisfied here, as are (7.5) and (8.9) with H(x) = ^Y17=i^ S x i 
and k(x) = (1 — C 2 )(^i H 1" x n ) X^j=i(l/(^ x j))- This function k(-) is non- 
negative, since we have assumed < (" < 1, whereas g(x) = (x\ + ••• + 
x n )(xi ■ ■ ■ x n )~ K . In particular, with £ = 1 we get 



(12.3) U(T, x) = Xl '-- Xn E 1 

Xl + --- + x n 



Xl (T) + ■ ■ ■ + X n (T) 
X\{T) ■ ■ -X n (T) 



[see Goia (2009) and Pal (2009) for a computation of the joint density of 
X\ (T) , . . . , X n (T) which leads then to an explicit computation of U (T, x) in 
(12.3) above, and shows that this function is indeed of class C > 2 ]. 

With C = 1 one computes Z(t) = U]=i( x j(°)/ X j( t ))^ therefore A(t) = 
(A(t)/A(0))"- 1 n"=i(^(i)M-(0)) as well as A t (t) = (X (t) / X (O))^ 1 ■ 
Ylj^ii^jd)/ 'a*j(0)) for i = 1, ... ,n. Both representations in (7.4) hold for 
the first hitting time of (7.1) in this case; whereas S = T = mini<j< n 7i as 
in (7.1H7.6), since L(t) = (1/Z(t)) = (X(t)/X(0)r U]=M (*)/^(0))- 

Both (9.12) and (9.13) hold for the example of (12.1) with h = n - 1, the 
first as equality; from the corollary to Proposition 2 and Remark 3, (9.16) 
holds. We recover the result of Banner and Fernholz (2008) on the existence 
of arbitrage relative to market of (12.1) over arbitrary time-horizons. 
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The diffusion process 2)(-) of (9.20) takes now the form 

(12.4) dYi(t) = Yi(t) dt + yjY$!)i?Jt) + • • • + Y n (t)) dWi(t). 

The conditions of Bass and Perkins (2003) are satisfied again, though one 
should compare the "weak drift" b.;(x) = x% > in (12.4), which vanishes for 
Xi = 0, with the "strong drift" bj(x) = k(x\ + • • • + x n ) for the the original 
diffusion X(-) in (12.2), which is strictly positive on [0,oo) n \ {0}. 

The corresponding Fichera drifts in (9.25), (8.8) are given by 2fj(x) = 

Xi - (xi H + x n ), 2fj(x) = C(x 1 H + x n ) - x i: and fj(x) > >fi(x) 

hold on {xi = 0} H {X^j^i^j > 0}> from Remark 2 we verify again that the 
diffusion £(•) of (12.2) takes values in (0,oo) n . 

In contrast, the new diffusion 2)(-) of (12.4) lives in [0,oo) n \ {0}, and 
hits the boundary O n of this nonnegative orthant with positive probability 
Q X [T < T] = 1 — U (T, x) for every T £ (0, oo). The positive P x -supermartingale 
Z(-)X(-) is a P x -potential, for every xG (0,oo) n . In this case, the three 
inequalities of (5.5) hold for every T £ (0,oo): the local martingales Z(-), 
Z(-)X(-) and Z(-)Xi(-), i = 1, . . . , n are all strict. 

The model (12.1) can be cast in the form (10.1) for the relative market 
weights, as a multivariate Jacobi diffusion process with dynamics d(ii(t) = 
(1 + C)(l - nm(t)) dt + y^t) dWi{t) - ia(t) Ek=i V^k(t)dW k (t), or 

(12.5) dm(t) = (1 + C)(l - nin(t)) dt + ^iM(t)(l-ia(t))dW}(t) 

with appropriate Brownian motions W|("), . . . , Wj(-). Thus, each compo- 
nent Hi(-) is also a diffusion on the unit interval (0,1) with local drift 
(1 + 0(1 — n u) an d local variance y(l — y) of Wright -Fisher type. Goia 
(2009) studies in detail this multivariate diffusion //(•) based on an exten- 
sion of the Warren and Yor (1999), Gourieroux and Jasiak (2006) study of 
skew-products involving Bessel and Jacobi processes. 

From (12.4), Y(-) := Yi(-) + ■ • ■ + Y n (-) satisfie s the stoch astic equation 
dY{t) = Y(t)[dt + dB(t)}, where B(-) := £™ =1 f Q y/Yj(t)/Y(t) dWj(t) is Brow- 
nian motion; thus Y(-) a geometric Brownian motion with drift, under Q x . 
The process v(-) = ■ ■ ■ ,u n (-))' of (10.2) is related to the auxiliary dif- 

fusion 2)(-) of (12.4) via Vi{-)=Yi{-)/Y{-). 

The dynamics of these fi(-)'s are easy to describe in the manner of (9.26), 
namely, dui(t) = sj Vi(t) dWi(t) — Vi(t) Ylk=i V u k{t) dWk{t), or in the nota- 
tion of (12.5): dui(t) = y/vi(t)(l-i/i(t))dWf(t). Then the Feller test [e.g., 
Karatzas and Shreve (1991), pages 348-350] ensures that each Ui(-) hits 
one of the endpoints of (0, 1) in finite expected time. Thus, all but one 
of the Yi(-)'s eventually get absorbed at zero; from that time X* [with 
E^ X (T*) < oo] onward, the only surviving nonzero component Y(-) behaves 
like geometric Brownian motion with drift; in particular, 2)(-) never hits the 
origin. 
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13. Some open questions. What conditions, if any, on the Markovian 
covariance structure of Section 8 will guarantee that 7?(-) of (11.1) never 
borrows from the money-market, that is, Y27=i x iDiU (T, x) < 0? That it is a 
portfolio, i.e., that ^?=i x iDiU(T,x) = holds? (See Remark 6 for a partial 
answer.) Or better, that tt(-) of (11.1) is a long-only portfolio, meaning that 
both this condition and A(G(T, x)e H(x) ) > hold? 

Can an iterative method be constructed which converges to the mini- 
mal solution of the parabolic differential inequality (1.1), U(0, ■) = 1 and is 
numerically implementable [possibly as in Ekstrom, Von Sydow and Tysk 
(2008)]? How about a Monte Carlo scheme that computes the quantity 
C7(T,x) of (9.22) by generating the paths of the diffusion process 2)(-), then 
simulating the probability Q X [T > T] that 2)(-) does not hit the boundary 
of the nonnegative orthant by time T, when started at 2)(0) = x € (0,oo) n ? 

How does U(T, x) behave as T — > oo? If it decreases to zero, then at what 
rate? 



14. Note added in proof. In the context of Proposition 1, and under 
the probability measure Q of Section 7.1, the processes Xi(-), . . . ,X n (-) are 
real-valued (do not explode) if and only if their sum X(-) as in (3.2) is real- 
valued. Now it is fairly straightforward to check from (6.6) that this sum 
satisfies the equation 

dX(t) = X(t)[d(M)(t) + dM(t)], 

where the continuous, Q-local martingale M(-) and its quadratic variation 
process (M)(-) are given, respectively, as 

M(t) : =E/ ( X^OO^OO J dW k (s), (M)(t) = ix'(s)a{s)n{s)ds. 

Thus by the Dambis-Dubins-Schwarz result [e.g., Karatzas and Shreve (1991), 
Theorem 3.4.6], for some real- valued Q-Brownian motion B(-) we have 



, < t < oo. 

u={M)(t) 



It is fairly clear form this representation that a sufficient condition for 
the total capitalization process X(-) to be real-valued, Q-a.e., is that this 
should hold for the quadratic variation process (M)(-): 

Q((M)(t)<oo,Vt€ [0,oo)) = l. 

In the volatility-stabilized model of Section 12.1 we have c%'(t) = 5ij/Hi(t) 
and thus (M)(t) = J27=i Jo ^(s) ^ s = so * n ^ s condition is clearly satisfied. 
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